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Abstract. A class of continuum models with a critical end point is considered whose Hamiltonian Tt[cj>,ip] 
involves two densities: a primary order-parameter field, <f>, and a secondary (noncritical) one, xp. Field- 
theoretic methods (renormalization group results in conjunction with functional methods) are used to give 
a systematic derivation of singularities occurring at critical end points. Specifically, the thermal singularity 
~ |t| 2 ~ a of the first-order line on which the disordered or ordered phase coexists with the noncritical 
spectator phase, and the coexistence singularity ~ |t| 1-a or ~ \t\@ of the secondary density (t/j) are derived. 
It is clarified how the renormalization group (RG) scenario found in position-space RG calculations, in 
which the critical end point and the critical line are mapped onto two separate fixed points Vqep- an d T^x 
translates into field theory. The critical RG eigenexponents of Vqep an d V\ are shown to match. 7-cep is 
demonstrated to have a discontinuity eigenperturbation (with eigenvalue y = d), tangent to the unstable 
trajectory that emanates from Pcep an d leads to V\. The nature and origin of this eigenperturbation as 
well as the role redundant operators play are elucidated. The results validate that the critical behavior at 
the end point is the same as on the critical line. 

PACS. 64.60.Fr Equilibrium properties near critical points, critical exponents - 05.70.Jk Critical point 
phenomena - 68.35.Rh Phase transitions and critical phenomena - 11. 10. Hi Renormalization group evo- 
lution of parameters 



1 Introduction 



Critical end points are widespread in nature. They oc- 
cur when a line of critical temperatures (or lambda line) 
T c (g), depending on a nonordering field g such as chem- 
ical potential or pressure, terminates at a line ga(T) 
of discontinuous phase transitions In the past 

decades plenty of experimental and theoretical studies 
have been made in wh ich crit ica l end p oi nts were e ncoun- 
tcred 




3 14,15^,16,17 



30T. Yet they have been rarely 
investigated for their own sake. This may be due to the 
physically appealing and widely accepted, though seldom 
carefully checked, expectation that the critical phenomena 
at such an end point should not differ in any significant 
way from critical phenomena along the critical line T c (g) 

§• 

However, Fisher and collaborators have pointed out re- 
cently §M|ll],|| that even the bulk thermodynamics of a 
critical end point should display new critical singularities, 
not observable on the critical line. For concreteness, let us 
consider the simple critical end point situation depicted 
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in Fig. [l| of a binary fluid mixture, whose critical line, 
is restricted to values g > g e and ends at the end point 
temperature T c = T c (g e ). We assume that the critical heat 
singularity at constant g > g e on the critical line can be 
written as A±(g) \T - T c (g)\~ a , where a > in d < 4 
dimensions. Fisher's assertion then is that the first-order 
phase boundary should vary as 



9*(T) « g^(T) 



(2-a)(l -a) 



(1) 



as t = (T - T c )/T c -> ±0. Here ^° S (T) is a regular back- 
ground term, and / should be equal to the universal 
(and hence g independent) ratio A + /A_. 

Equation (|l|) has been derived by using general ther- 
modynamic arguments in conjunction with the phc- 
nomenological theory of scaling A straightfor- 

ward extension of such reasoning reveals that the ther- 
modynamic density conjugate to g, which in the present 



1 For the sake of simplicity, we shall henceforth refer to the 
critical line A briefly as the 'A-line', a term usually reserved for 
the critical line of normalfluid-to-superfluid transitions of He. 
Since we shall explicitly consider only the case of a scalar order 
parameter cj> (even though parts of our analysis can be general- 
ized in a straightforward fashion to more general cases with an 
multi-component order parameter), no confusion should arise 
from this slight abuse of terminology. 
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Fig. 1. Upper drawing: Schematic phase diagram showing, 
in a field space, portions of the coexistence surfaces p and a 
for a nonsymmetric binary mixture with a liquid-liquid critical 
line A and a critical end point CEP (full circle). The lower 
drawing represents the projection of the phase diagram on the 
gT-plane. Along the triple line the ordered phases f3 (A-rich) 
and 7 (A-poor) coexist with the spectator phase a (vapor). 



case can be identified as the total particle density gtoti 
displays nonanalytic behavior of the form 



sing 
ftot 



ul 



\tf + V£ |i| J 



as the end point is approached along the liquid side of the 
liquid- vapor coexistence curve from t > or t < |T^,[h|. 
Here g r t °f (T) is regular at T c . Further, the amplitude U± 
vanishes for symmetric binary fluids whose properties are 
invariant with regard to simultaneous interchange of the 
two constituents A and B and their respective chemical 
potentials fiA and /ib- 

The |£| 2_a singularity of (jl]) has been confirmed for ex- 
actly solvable spherical models pi and checked by Monte 
Carlo calculations the W~° singularity of (§) is 

consistent with the jump in the slope of Ptnt(T) found in 
mean field and density functional calculations ||19|,Pl| and 
has also been seen in Monte Carlo simulations ll 3L |1 4, . 



Despite these advances the present state of the the- 
ory is anything but satisfactory. As can be seen from the 
thorough discussion given in Ref. H, a number of chal- 
lenging problems exist. First, a systematic derivation of 
the above-mentioned singularities within the framework 
of the modern field-theoretic renormalization group (RG) 
approach is lacking. We are aware of only a single field- 
theoretic analysis of critical end-point behavior that goes 
beyond Landau theory: the (e = 4 — rf)-cxpansion study 
of a scalar model with a negative 6 term by Ziman et 
al. H . In their one- loop calculation they found that both 
the critical line and the end point were controlled by the 
same, standard 0(e) fixed point. However, this model has 
very special properties: The first-order phase boundary 
does not extend into the disordered phase; as the critical 
end point is approached from the disordered phase, the 
order parameter (</>) becomes critical and exhibits a jump 
to a nonvanishing value upon entering the ordered phase; 
no critical fluctuations occur in the ordered phase. Thus 
the model clearly does not reflect the typical critical end 
point situation in which the two-phase coexistence surface 
p bounded by the critical line T c (g) meets the spectator 
phase boundary <x in a triple line (see Fig. 

Second, the above RG scenario for critical end points 
clearly differs from the one encountered in position-space 
RG analyses of lattice models ^,E5|. These yielded the 
RG flow pattern depicted in Fig. in which the crit- 
ical line is mapped onto the fixed point while the 
critical end point is described by a separate fixed point 
■pCEP- Since separate critical fixed points normally rep- 
resent distinct universality classes of critical behavior, it 
would be natural to expect that the corresponding criti- 
cal RG eigenexponents yt = \ jv and ?//, = Ajv take on 
different values at these fixed points. In Ref. p2| , |25| ] these 
values were found to match. Further, Vq EP was found to 
involve an additional relevant eigenvector with eigenexpo- 
ncnt y = d, characteristic of a discontinuity fixed point 
& 




CEP 



Fig. 2. Schematic RG flow pattern found for critical end 
points in position-space RG calculations J2^,^5|. The fixed 
point "Pcepi representing the critical end point, has two rele- 
vant eigenperturbations whose eigenexponents match those of 
the fixed point V\ , and a discontinuity eigenperturbation along 
the trajectory connecting the two fixed points, with eigenex- 
ponent y — d. 
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Whether this RG scenario — two separate fixed points 
V\ and "Pqep with matching critical spectra — is correct, 
remains to be seen. According to Fisher and Barbosa S, 
it need not be the invariable rule, even though it should 
normally be expected. If the latter is true, then one ought 
to be able to corroborate this RG scenario by means of a 
systematic derivation based on the field-theoretic RG ap- 
proach. This should also yield a discontinuity eigenexpo- 
nent at "PJep , which the field-theoretic analysis of Ziman 
et al. || was unable to give. 

A natural way to tackle this problem is to consider 
models involving two fluctuating densities, namely, a (pri- 
mary) order parameter field (p> and a secondary (noncriti- 
cal) density ip. The reason should be clear: Across p, the 
order parameter has a jump singularity, but across er, an- 
other density — the thermodynamic density conjugate to 
g — is discontinuous. Obviously, a proper thermodynamic 
description requires two thermodynamic densities, in ad- 
dition to the thermodynamic field T. Models of this sort 
have been investi gate d by means of Landau and density 
functional theory |l9|,|(| and by Monte Carlo cal- 
culations C5). Unless long-range interactions have been 
included also, they may be viewed as continuum versions 
of the Blume-Emery-Griffiths (BEG) model [||. 

The continuum models with short-range interactions 
we shall consider here arc of this kind but more general 
than the previously studied ones. They have been chosen 
in such a way that they can also be used to investigate in- 
terfacial critical phenomena such as the critical adsorption 
of one of the two components of the fluid at the interface 
between the liquid /3~f phase and the spectator phase a 

@,HH,0. 

In the present paper^J we shall confine our attention to 
bulk properties, leaving the study of critical adsorption, 
with the challenging question of whether the spectator 
phase a may be replaced by a hard wall (the 'wall as- 
sumption' of Ref. [till , implicit in virtually all work on 
critical adsorption at interfaces), to a forthcoming paper 

Our principal aims are the following. First, using field- 
theoretic means, we wish to give a systematic derivation of 
the coexistence singularities (Q) and (Q) beyond the level of 
Landau theory. Our second goal is to show that the critical 
behavior at the critical end point is the same as along the 
critical line A, inasmuch as universal properties are con- 
cerned. To this end we shall have to clarify whether the 
RG flow pattern found in position-space RG calculation 
p2|,p5| prevails in our field-theoretic analysis. In particu- 
lar, we shall verify the appearance of a g-dependent scaling 
field with RG eigenexponent y = d, and explain its origin 
and significance (cf. Ref. pg|]). 

An outline of the paper is as follows. In the next sec- 
tion, we introduce appropriate continuum models both for 
the case of symmetric and of nonsymmetric binary flu- 
ids; we justify their choice on the basis of phenomenolog- 
ical arguments, and show that they can be obtained from 
the BEG model by means of a Kac-Hubbard-Stratonovich 



transformation and subsequent continuum approximation. 
In Section |^, we utilize the Landau approximation to an- 
alyze these models; we verify that for suitable choices of 
their parameters, phase diagrams with the correct qualita- 
tive features result, determine the Landau-theory analogs 
of the coexistence singularities (|l|) and (|J), and work out 
a number of perturbative results that will be needed in 
later sections. The field-theoretical analysis beyond Lan- 
dau theory is taken up in Section ^[ Section || summarizes 
our main results and conclusions. Finally, the Appendix 
contains some details of the mapping of the BEG model 
to the employed continuum models. 



2 Models 

2.1 Phenomenological considerations 

The models we wish to consider have a Hamiltonian of the 
form 

W[0,V]=Hi[0]+H 2 [V']+Hi 2 [0,V] , (3) 

where cp{x) and ip{x) are scalar fields on a bounded region 
fl of the (i-dimensional Euclidean space M. d . It is under- 
stood that the limit fl | R d is taken; since we shall only be 
concerned with bulk properties in the present paper, one 
can imagine fl to be a ci-dimensional hypercube with pe- 
riodic boundary conditions whose linear extension is infi- 
nite. While 4>, the primary order parameter field, becomes 
critical, the secondary field, ip, is noncriticalj^j 

To simplify matters, we shall ignore long-range interac- 
tions. In real binary fluids, long-range interactions of the 
van der Waals' type are normally present. While these 
play a crucial role for — and would have to be included in 
a study of - 



■ wetting phenomena |4l) , they can be trusted 
to be irrelevant in the RG sense for critical phenomena. 
Anticipating that a gradient expansion be made, we pre- 
sume that the Hamiltonian is a local functional depending 
on <p, ip, and their gradients. 

Consider first the case of a hypothetical binary fluid 
whose properties are symmetric with regard to inter- 
change of the two constituents <r = A, B. If the interaction 
between ^ and particles is described by a pair potential 
[/„/(?-), then this symmetry is realized if Uaa = Ubb and 
f-A = Mb, where /^a,b are the respective chemical poten- 
tials. For such symmetric binary fluids, one can identify (p 
and ip (up to convenient normalization factors) with the 
difference qa — £b and sum qa — Qb of the total local 
mass densities £>a,b, respectively. The Hamiltonian must 
be even in </>; that is, TL, Ti\, and TL12 must satisfy the 
condition 

H[-^]=H[4>,ip]- ( 4 ) 



2 A brief report of parts of the results presented here has 
been given in Ref. |p8|; see also Ref. |39||. 



3 The term 'noncritical' is not meant to imply that the cor- 
relation functions of the ip field do not display any critical 
singularities. Owing to the coupling between between ip and 
4>, correlation functions of the ip field also display critical be- 
havior, as should be clear and will be shown explicitly below. 
Only if this coupling vanishes, are the correlation functions of 
ip those of a massive field theory. 
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Table 1. Canonical momentum dimensions for d — 4 — e 



<P 




"2 


<7 4 


ft 


^ B 


1-1 


2 — - 

Z 2 


2 


C 


S — - 

J 2 


-2 










d 21 


en e 2 i,/2i 





e-4 


9 — 1 
z 2 


1 


c: 
2 


-1 1-2 



On the other hand, terms even and odd in ip are allowed. 

We shall keep in 7ii monomials up to fourth order in 
<j), in Ti.2 those to the same order in ip, and in 7ii 2 terms up 
to order cp 2 ip. Further, contributions of higher than second 
order in the gradient operator V will be discarded. It is 
convenient to make an cc-indepcndent shift 



tp{x) -> tp(x) , 



(5) 



choosing 'P such that the coefficient of the ip 3 term in 
TL 2 vanishes. The requirement of invariance under the Eu- 
clidean group E(d) gives further restrictions, forbidding, 
in particular, terms linear in V such as J n ((pVip ~ ip'Vcf)- 
Noting that J^ipAcp 2 equals j Q (p 2 Aip up to boundary 



in 

terms, we arrive at 



Hx[ct>] = [ 
Jn 



.Q 



(6) 



(7) 



and 



Hv»[<t>M = f 
Jn 



dn (j> ip + 4> 2 i> + en <p ^ 



with 



h = dn = en = 



(8) 



(9) 

At the Gaussian fixed point at which <p is critical while ip 
is noncritical, the coefficient a 2 vanishes, but b 2 remains 
nonzero. Since the coefficients A and b 2 could be trans- 
formed into unity by a change of normalization of <p and 
ip, it is appropriate to take A and b 2 as dimensionless |j42f . 
The momentum dimensions of the other parameters are 
listed in Table 0. 

If the restriction = Mb is lifted by turning on a 
chemical potential difference = — Mb while the AB- 
symmetry of the microscopic interactions persists, then 
property (0) will be lost, the interaction constants in (0) 
will no longer vanish, and a <p 3 term is to be expected in 
Til. The latter can be transformed away by making a shift 



(*) 



(x)+<P 



(10) 



analogous to (||) . The coefficients of all contributions that 
are odd in <p must change sign as /i_ is reversed. This prop- 
erty carries over to <P. Hence these interaction constants 
as well as # vary ~ /z_ for small 



In the general case, in which the interactions are not 
AB-symmctric, the interaction constants of the terms of 
TL that break the cf> — ► — <f> symmetry are no longer odd in 
/i_ and hence do not vanish for //_ = 0. The same applies 
to the shift <P. However, the Hamiltonian specified by (H) 
and remains appropriate. The main difference is 

that the coefficients dn and en do not vanish and we 
must consider nonzero values of the ordering field h. 

Our choice of monomials retained in 7i[4>,ip] requires 
some comments. The ^-dependent part "Hi is the standard 
<p A Hamiltonian, comprising all those monomials whose 
coefficients have nonnegative momentum dimensions for 
small e = 4-d > (cf. Table |) and are not redundant (as 
is the <p 3 term) |l3| . If besides Hi only the terms quadratic 
and linear in TL2 and the one cx d 2 \ in TL\ 2 were kept, then 
the Hamiltonian H. would agree with the one utilized in 
the definition of the familiar dynamic model C |E3 . 

The ip 4 term with 64 > has been included for two 
reasons. First, if b 2 < 0, it is needed for stability. Second, 
we ought to be able to obtain a spectator phase boundary 
a from our model. Therefore it is natural to require the 
model to yield such a coexistence surface already in the ab- 
sence of coupling to the primary field <p (i- e., for TL\ 2 = 0). 
The obvious analog of a for this case (and given 64 > 0) is 
the plane g — 0; within Landau theory, the two (spectator 
and 'liquid') bulk phases coexisting there for 62 < cor- 
respond to the expectation values (ip) = ± \Z\b 2 \/b4. Note 
that, rewritten in terms of the shifted and rcscalcd field 



ipd 



H 2 takes the form 




(11) 



Thus its V± 



■9 fa 



(12) 



term is positive and independent of 64. The 
coefficients of the ip± and ip± terms in ( |l2| ) (which have 
negative momentum dimensions) vanish as 64 — > 0. If we 
let these coefficients as well as the one oc B of the (Vip±) 2 
monomial (whose momentum dimension is also negative) 
approach zero, and set the field g to the value g = ±0, 
then the rp± -dependent part of TL 2 turns into precisely 
this quadratic part, namely the Gaussian Hamiltonian 



n G [ip ± ] 



n 



fa\ ■ 



(13) 



For vanishing Tli 2l this is a fixed-point Hamiltonian in 
the space of "0-dependent Hamiltonians. We refrain from 
a more detailed discussion of RG issues here, reserving it 
to later sections. Note, however, that when the previously 
mentioned gradient, cubic, and quartic terms of H. 2 are 
dropped together with all contributions in Hi 2 other than 
the one oc d 2 i, then the total Hamiltonian TL reduces to 
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the one of dynamic model C. A trivial, but useful, obser- 
vation is that the constant part of the Hamiltonian ( p"2|) 
has a g- dependent part. Any constant in the integrand 
of a bulk Hamiltonian couples to \f2\ = J^, the volume, 
and hence scales like a relevant scaling field with scaling 
index d under RG transformations. In Wegner's terminol- 
ogy jl3|], such a scaling field is called special rather than 
relevant, because of its well-known relationship with the 
analytic part of the bulk free energy density. In any case, 
we see already at this stage (where the coupling to the or- 
der parameter field <fi has not yet been taken into account) 
a scaling field emerging that is odd in g and has scaling 
index d. 

As stressed earlier, 64 must not simply be set to zero: 
It must have a positive value in order that our model may 
yield a bulk phase diagram with the correct topology, for 
appropriate values of the other parameters (cf. Sec. ||). 
The coefficient B, on the other hand, can be set to zero in 
essentially all of our subsequent analysis. Our main reason 
for including this term is our intention of employing the 
model in a forthcoming paper (4(J in a study of critical 
adsorption of binary fluid mixtures at their interface with 
the spectator gas phase. The coefficient B serves to pro- 
vide a scale (correlation length) on which kink solutions 
ip |_(z) of the Ginzburg-Landau equation for TI2, connect- 
ing the asymptotic bulk solutions =F \/ 1 62 1 / ^4 at z = ±00, 
vary along the direction perpendicular to the interface. 
As B — ► 0, these kinks become sharp steps. One does not 
need to be concerned about the vanishing of this scale here 
since we shall only deal with bulk properties below. 

Likewise, the Laplacian terms cx en and cx e2i as well 
as the term oc /21 (all of which have negative momentum 
dimensions) have been introduced primarily with a view 
to our intended analysis of interfacial problems , and 
can be ignored in the sequel. However, the coupling term 
cx g?ii, which is also not present in the Hamiltonian of 
dynamic model C, must not be ignored in the general 
nonsymmetric case. As the reader might anticipate and 
will see below, it plays an important role in producing the 
\t\P part of the coexistence singularity (0). 



2.2 Derivation from the Blume-Emery-Griffiths model 

A familiar lattice model for binary mixtures that exhibits 
a critical end point is the BEG model [^3[Q|. This is a 
classical spin-1 lattice model with the Hamiltonian 

Hbeg [S] = - J2 [J Si Sj + K Si Sj 

+ L(StSj + SiS?)] 
-J2( HS i + DS i)> Si = 0, ±1,(14) 



where indicates summation over nearest-neighbor 

pairs of lattice sites. Initially it was proposed to simu- 
late 3 He- 4 He mixtures; in this interpretation, a 3 He atom 
at site i corresponds to Si = and a 4 He atom to Si = ±1 



[ p3[ , so there is one and only one helium atom at each lat- 
tice site, and the model makes no allowance for vacancies. 

However, the model may also be interpreted in a 
distinct — and for our purposes more appealing — fashion 
as a lattice-gas model for a binary fluid ^5|. Then i is 
a label for microscopic cells that can hold at most a sin- 
gle A or B particle (atom or molecule), and the states 
Si = +1, —1, and correspond, respectively, to the cases 
'cell i is occupied by an A particle', 'cell i is occupied by 
a B particle', and 'cell i is empty'. Thus the variables H 
and D represent odd and even linear combinations of the 
chemical potentials ma and Mb p5|: 



and 



H 



D 



MA - MB 

2fc B T 



MA + Mb 
2k B T 



(15) 



(16) 



Without loss of generality, we can take L > Jlj| . We also 
assume that J > and K > 0. Thus, as J increases, the 
tendency for phase separation in the liquid state grows, 
and larger values of K correspond to a stronger drive of 
condensation of the mixture. 

On sufficiently large length scales, the physics of the 
BEG model should be described by the continuum models 
introduced in the previous subsection. Hence one ought to 
be able to derive the latter from it. This may be achieved 
in much the same way as one can map the Ising model 
onto a continuum field theory, utilizing a Kac-Hubbard- 
Stratonovich transformation (see, e.g., Ref. [Q). Details 
are given in the Appendix, where we show that the grand 
partition function of the BEG model, 



-BEG 



E ■ 

{5i=0,±l} 



-Wbeg [S] 



(17) 



can we rewritten exactly as an integral over continuous 
fields 4>i and ipii that is to say, as the partition function 
of a lattice field theory. One has 



-BEG 



-h{J,L,K) 



n 



di/ii 



-Hm\<t>,i>] 



(18) 

Th e exp licit form of this Hamiltonian TC\[ t [<fi,i(>] is given 
in (183) of the Appendix. As is shown there, making a 



continuum approximation then yields a Hamiltonian of 
the form specified by (0) and (0)-(||). 



3 Landau theory 

Landau theory is of value not only because of its simplic- 
ity and ability to reproduce essential topological features 
of the phase diagram, but also because it serves as the 
starting point of studies based on RG improved perturba- 
tion theory. Our main goal here is to convince ourselves 
that our model indeed yields a bulk phase diagram with 
a critical end point and the topology illustrated in Fig. 
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providing the model parameters are in the appropriate 
range. 

Since we shall only consider translationally invari- 
ant states here, we can take <f> and ip to be position- 
independent in the rest of this section. Hence all terms 
of TL involving spatial derivatives do not contribute. If we 
drop these terms, then its part even in <f> becomes equiva- 
lent to the Hamiltonian investigated in Ref. jl9] via Lan- 
dau theory. Thus, for the symmetric case, the results may 
be partly inferred from this reference. 

In the Landau approximation, the grand potential (per 
volume and k^T) is given by 



with 



A(a, b, d, h, g) = inf V(<f), ip) , 

9>V> 



(19) 



(20) 



where a, b, and d stand for the sets of parameters {02, 04}, 
{62, 64}, and {dn, ^21}, respectively. From (0)-(§) we find 



V(0,V) 



(a 2 + d 21 i>) + ?± 4 - H 



|^ + ^4 + du ^_^_ (21) 



For the sake of thermodynamic stability, 04 and 64 are 
assumed to be positive. 

In place of g we shall occasionally use the conjugate 
density 

* = ( 22 ) 

°9 a ,b,d,h 

as independent thermodynamic variable, utilizing the 
mixed field-density representation (a, b, d, h, ip) instead of 
the field representation (a,b,d,h, g). The thermodynamic 
potential associated with the former is defined by 



B(a,b,d,h,$) = w£{V{<f>,$) + g$} 



(23) 



3.1 Symmetric case 

We first consider the symmetric case, setting dii = 0. 
Owing to the implied invariance with regard to (</>, h) — > 
(— <j), —h), we may restrict ourselves to values h > 0. Since 
a sign change of d 2 \ can be compensated by simultaneous 
sign changes of tp and <?, we may furthermore presume that 



(hi < 



(24) 



From the analysis presented in Ref. |l9| it is clear that a 
variety of different types of phase diagrams (with critical 
points, critical end points, tricritical points, or special tri- 
critical points) can be obtained, depending on the chosen 
range of the parameters a 2l b 2 , d 2 \, and g. Our aim here is 
not to explore all these possibilities; rather we shall focus 
on the case of critical end points, and choose the values of 
these parameters accordingly. 



The equilibrium densities must satisfy the (classical) 
equations of state 



0, 



(25) 



d(f> dip 
which in the present symmetric case become 

(a 2 +d 21 il))4> + a i 4? = h , (26) 



62 i> + &4 i> + 



3 , rf 2i a2 



(27) 



From the cj> — > — <f> symmetry (||) it is clear that the 
coexistence surface <x of the ordered (</> ^ 0) states must lie 
in the h — plane. Upon setting h = 0, the first equation 
of state, (0), can easily be solved for to determine that 
value ^min(V') & t which V((f>, tp) becomes minimal for the 
given value of ip. One finds 

r , if a 2 + d 21 iJj > 0. 

Ml>) - I r. - „ (28) 

The critical value of ip at which the bifurcation occurs, 

4>X = , (29) 

is the solution to 



2 



V(0,V) = 0, 



(30) 



a condition that must be fulfilled whenever <fi is critical at 
= and hence on the A-line. 

Upon inserting ( ^8| ) into the right-hand side of (p3|), 
we obtain 

^2 ,2 . ^4 ,4 



B(a,b,c,0,1>) = + 



a 2 



ipx I 4a 4 



*21 



(iP~^x) 2 ,(31) 



where 9(.) is the step function. 

Since we wish the density ip to be positive on the A-line 
we take 

a 2 > ; (32) 

together with the inequality ( |24| ) this ensures that ip\ > 0. 
Under these conditions the trivial solution 4> m - ln = holds 
in the regime ip < ip\. Thus the equilibrium densities ip a 
and ippy of the disordered vapor and liquid phases are 
solutions to the = analog of (p7|): 



6 2 ip + 64 ip 3 = g 



(33) 



If b 2 > 0, there exists a unique real solution, which is 
positive or negative, depending on whether g > or g < 0, 
namely 



1 



Vo(b 2 ,b 4 ,g) = - (r ) 



V3 _ £ (ro) -l/3 



(34) 
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with 




(35) 



It has the power-series expansion 

" big 3 , 3blg 5 



Voifo > 0,64,5) = ^ ~ & | 



+ ^ + 0(.9 7 ), (36) 



where the subscript is to remind us that the solution 
vanishes for 5 — 0. For 62 < 0, ( p3[ ) has three real solutions 
if 



\9\ <9< = 



_ 2\b: 



13/2 



(37) 



3V36I 

and a single real one if I5I > <?<• We denote the ones that 
turn in the limit g —> into the nontrivial g 7^ solutions 
=t VI ^2 1/^4 as 77-1- (62 , 64, 3) • They have the property 



V-( b 2,b 4 ,g) = -774.(62,64,-5) 
and the power-series expansion 



(38) 



fe(6 2 < 0,64,5) = ±\r-r 

04 




64 5 3 _ 105 6^ V 



2|6 2 | 4 1 128|6 2 |"/2 
3 6 2 n 5 

S + 0(5 6 ). 



(39) 



Further, 77+ is given by the analog of the result (34) for 770 
one obtains through the replacement 62 — * — |62]7both in 
( |34| ) and the expression ([35]) for r . 

In terms of these solutions, the densities of the disor- 
dered fluid state become 

= 0, ip = 770(62, 64, 5) , for 6 2 > , h = , V < i>X , 

(40) 

while the densities of the disordered vapor phase a and 
the disordered liquid phase /3"f read 

(0 a ,Va) =(0,^(62,64,5)) . \ h u 2< ^ 

for < h = , (41) 
Wti^) = (°> 7 7+( & 2,6 4 ,5)) I V < V'a • 

In the latter two cases, the solutions 77+ (/?7 phase) 
and ?7_ (a phase) are the thermodynamic stable ones if 
5 > and g < 0, respectively; the remaining two real 
roots of ( |34| ) one has if I5I < g < correspond to metastable 
and unstable states. 

At 5 = 0, these states coexist down to that value of 62 
at which ipp^ intersects the line tp = ipx> which is 



-64^ 



d.2i 



provided 



-d 2 i 



2 \J 04 64 



(42) 
(43) 



As we shall show below, this condition guarantees that no 
tricritical point appears. The resulting phase diagram is 
displayed in Fig. |3j in a mixed field/density representation, 
and in Fig. | in a field representation; it corresponds to 
case (e) of Ref . |^9| .^| There is a liquid- vapor critical point 
at 62 = ip = and a critical end point at 62 = 62 e , ip = ipx- 
The boundaries of the a-0'y coexistence region for h 2e < 
6 2 < are given by the 5 = solutions r]± of (|39|). 




Fig. 3. Phase diagram of the symmetric model in the Landau 
approximation. The diagram is shown in a mixed field/density 
representation and corresponds to the following choice of pa- 
rameter values: d2 = 5, 04 = 64 = 1/6, d,2i = —1/2, dii = 0, 
giving ipx = 10 and &2e = 50/3. 



In order to determine the boundaries of the region of 
coexistence of the a, /3, and 7 phases, one generally must 
resort to numerical methods. However, for 62 < 62c the 
value 5 takes on the triple line, g Gl is small, as is the order 
parameter <fi for ip > ip\ in the vicinity of the critical end 
point. This enables us to determine the phase boundaries 
in this regime in a perturbative manner. Note, first, that 
the equilibrium densities ipp and ip-y (= ipp for h = 0) can 
be written in terms of 



<MflO = 77+02,64,5) 



and expanded as 



(44) 



(45) 



with 
and 



#, 7 = v+ H>2, 64, 5 - - ^21 4>p n 

= 1> + (g)-fl>+'(jg)^4„ 

+ ^ + "(5)^^ 7 + 0(^ 7 ), (46) 

1>+'(g) = [364^(5) -N]" 1 (47) 

6 64 ^+(5) 



</>+"(<?) = - 



364^(5) -N] 



3 • 



(48) 



4 In order to check the consistency with Ref. [|19[ , one should 
identify the variables q, a, /j and A utilized there with — ip, 02, 
—g, and 62, respectively, and set an = g^i = 62 = 64 = 1. 
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Fig. 4. The analog of Fig. ^ in an a-igh field representation, as 
obtained by solving the equations of Landau theory through 
numerical means. For 62 we have chosen the value 62 = —6; 
the remaining parameters have the same values as in Fig. 
i.e., 04 = 64 = 1/6, d,2\ = —1/2, and tin = 0. Only the neigh- 
borhood of the critical end point is shown. 



Inserting the expansion (46) into the equation of state (|26j 
with h = gives 



[a 2 + d 21 ip+{g)] - 



d 2 

04 - —ijj + '(g) 

4x 



+ ^V+"(ff)^ 5 -O(0 7 ) • (49) 
We seek nontrivial solutions of this equation for 
a 2 [ip+{g)} = a 2 + d 21 ip+(g) < . 



If the coefficient of the 6 term in 



is positive, 



a 4 [^ + ( 5 )] = a 4 -^ + '(5) >0, 



(50) 



(51) 



then the equation ( |49| ) for (j> can be solved in a straight- 
forward manner to obtain 



V 04 ^+ 3) 



giving 



d2i \a 2 [ip+(g)}\ 
2 a 4 [V+(5)] 



(52) 



(53) 



The above analysis suggests that 04 [ip+(g)] plays the 
role of an effective 4> A interaction constant. To verify this 



one can calculate d 4 V[0, p ^ ((/))] /defy 4 \^ =Q along a path 
ipf3~{ (4>), where ^g i7 (0) is given by (p5|), with the replace- 



ment 



A straightforward calculation gives 



d 4 V 



6 CL4, [tp+(g)] 



(54) 



6=0 



where the equation of state (|26|) and the results ( |47|) a nd 
( p8| ) for the derivatives of V>+ (ff) were used. Hence (|5l|) is 
the usual condition for the absence of a tricritical point at 
= 0. At the critical end point it becomes 



CI4 > 



4I6 



2c 



d 2 i 



(55) 



which in turn implies our initially stated condition (43) if 
d 2 i < and tpx > 0, as assumed. 

The location of the triple line, g ai can now be deter- 
mined from the coexistence condition 

A a (a,b,d,0,g a ) = -4/3,7 (a, b > d > °>9a) (56) 

for the grand potentials 

A p (a,b,d,h,g) = V{(/) p ,ip p ) , p = a,j3,^. (57) 

Upon substituting the above results for (0/3, 7 , ipf},~{) to- 
gether with {4> a ,Tpa) from (41) into the latter equation, 
the former one can be solved beneath the critical end point 
to obtain 



[W)f 

8 di(ip) i/j 
a 2 + d 2 \ 



O 



([«)] 3 ) 



V>=V>+(0) 



8 (a 4 



4|fc 2 | 



O(dl) 



(58) 



(59) 



< 0. 



l 2.± 



\t\ near the crit- 



for a 2 [\/\b^\fb^\ = a 2 + d 2 i 

The latter quantity varies as a 2 
ical end point, while the denominator in (|38|) approaches 
a nonvanishing constant there. Remembering that g a = 
for a 2 > (i.e., £ > 0), we arrive at 



dt 2 







t=o± 



4a 4e ^/|62e|/&4e 



for t = 0+, 
for t = 0", 



(60) 



where the subscript e denotes values at the critical end 
point. Hence d 2 g a /dt 2 has indeed a jump singularity at 
the end point, as found in other analyses of the mean-field 
type |lj|^0],[l^,|3l[] and is in conformity with the a = 
analog of the predicted singularity (|l|) . 
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Likewise, the singularity^] of the liquid (/?7 or /3, 7) 
branch of the a-(3j and a-fl-j coexistence regions, 



^21,e |«2 



07/(9,7 



416 



2e| «4e 



\t\6(-t) 



(61) 



resulting from the second term in ( |53| ) is the mean-field 
analog of the one ~ |^| 1_Q m (|^); it means that in Landau 
theory the first temperature derivative of ipn is discon- 
tinuous at the critical end point (cf. Refs. [^[ij). 

The curve in which this portion of the coexistence 
boundary intersects the 62V' plane — i.e., the right branch 
of the liquid phase boundary in Fig. ||, described by the 
mapping 62 = > ipp-y and 62 1— * ipp.-y f°r > 62 > b 2e and 
02e > 02, respectively, — also has a discontinuous slope at 
the end point. One finds 



<%2 



-\d 2 



b 2 =b 2e +0 



2 0,2 64 



(62) 



and 



db 2 



-\d 



21 



l 21\ 



20264 2a2 64 (4a| 04 64 — |c?2i| 4 ) 
2a 2 04 \d 2 i\ 



4 a 2 . 04 64 — |d2i 



(63) 



3.2 Nonsymmetric case 

We now turn to the case of a nonsymmetric critical end 
point, assuming that d\\ > 0. Owing to the implied ex- 
plicit breaking of the (p — > — (p symmetry, the coexistence 
surface p should neither be located in the h = plane 
nor be parallel to it. Hence we must consider nonvanish- 
ing values of h from the outset. Instead of ( |26| ) and (p7|), 
the classical equations of state (Ea) now read 



(64) 



(a 2 + d 21 ip) <p + dn ip + a 4 (j) 3 



6 2 V + ^4^ 3 + dii(/)+%0 2 =g. (65) 



and 



Using the solutions 770, ± introduced in (|3^ 
l), the former can be solved for <fr to obtain 

(f) = rjo,±( a 2 + d 2i4>,a4 1 h - dnip) , 



(66) 



where the subscript is to be taken whenever 02 + d 2 \ ip > 
0, while the choices ± apply to the cases 02 + d 2 \ ip < 

equation of state, (|64j), for ip gives 



with h—dn ip ^ 0, resp ectively. Likewise, solving the other 



i> = Vo,±[ 62,64,3-^11 



d 2 \ 
2 



(67) 



5 To obtain the leading temperature singularity given here, 
one may replace flo-(t) by its limiting value g e = g a (t — 0); the 
\t\ 2 singularity ( |60[ ) of g a (i) produces subleading nonanalytic 
contributions to ippy/p,y. 



Since we shall continue to take 62 < 0, the subscripts ± 
are the appropriate choices, depending on whether 



S(4) = g-d n (f>-(dai/2) 4? 



(68) 



is positive or negative. 

From the above equations one can read off that the so- 
lutions </> m in(V0 given in ( p8|) now apply on the hyperplane 
h = dn ip. If 02 + d 2 \ ip > 0, then any thermodynamically 
stable state on this hyperplane must have = 0. Fur- 
thermore, the liquid-gas critical point is seen to remain at 
b = g = h = 0. 

Since <p\, the equilibrium value <p takes on the A-line, 
generically does not vanish, the fields p and ip do no longer 
decouple there. To see this, let us compute the Hessian of 
V, 



V (2) (0,^) 



o 2 v 



d<p d<p d<f> dip 

d 2 V d 2 V 

d<p dip dtp dip 



(69) 



at a solution (p = (p c \, ip = ip c \ of the equations of state ( |64| ) 
and (|65| ) ("classical solution"). This quantity is nothing 
but the Fourier transform 



r (2) (,) = 



fd d x 12 r flv (x 1) x 2 ;)e i "-^ 



(70) 



at momentum q = of the familiar Ornstein-Zernicke 
expression for the two-point vertex function 



r ilv (x 1 ,X2) 



s 2 r[p=p ch ip=ip cl ] 

5fj,(xi) 8v{x 2 ) 



(71) 



Here the indices [i, v take the values 0, ip. Further, x\ 2 
means the deplacement vector x\ — x 2 , and r[<p,ip] = 
H[4>,ip] in the Landau approximation used here. 
As is borne out by the result 



V (2) {ct> ch ip c 



a 2 + d 2 i ipd + 3 a 4 4> 2 C \ dn + d 2 \ p c \ 
du + d 2 i (p c i b 2 + 3 64 i>l\ 



(72) 

the Hessian now is generically nondiagonal, even in the 
= plane (where it is diagonal in the symmetric case) . 

A principal axis transformation with the orthogonal 
transformation matrix 



U 



cos "i? sin ?9 
— sin i? cos 1? 



(73) 



yields the diagonal matrix 

diag (Ai , A 2 ) = U T ■ V {2) (<p ch ipci) • U (74) 
with the eigenvalues 

A* = V00 \ V ^ T \ V flW - Vu? + 4 V^ a , (75) 
where the angle 1} is given by 

2 V^tf, 



tan 2$ 



(76) 
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We write the associated eigendensities as 



Pi 

P2 



<Pcl 



(77) 



On the A-line, the Hessian V^(0a, ip\) must have one 
vanishing eigenvalue, Ai, and a positive one, A2: 



Ai(0a,^a) = O, A 2 (0 a ,^a)>O 



(78) 



In order that V((f>, tp) > V{4>\, ip\) for small deviations <pi 
with <p2 = 0, the third derivative of V along the eigendi- 
rection ipi, 



d 3 V 

dpi 



= 6 0,4 (f> c \ cos 3 $ — 3 g?21 cos 2 i9 sin$ 



— 6 &4 ^ci sin 3 1? , 
must vanish on the A-line: 
8*V 

dpi, ' 
while the corresponding fourth derivative 

cl 



dpi 



6 04 cos + 6 64 sin 



(79) 



(80) 



(81) 



has to be positive there. The latter is guaranteed by our 
assumptions that 04 > and 64 > 0. 

Finally, the analog of the requirement (^lj) that the 
effective 4 coupling constant be positive becomes 



dpi 



V 



a 4 ( 



cl ,Vci)>0, (82) 



where the derivative on the left-hand side is along a path 
(pi h- > ^'"(Vi) through (4>ch"tpci) on which <^2 takes that 
value '/'^'"('/'l) which minimizes V for given ip\. Exploiting 
the fact that 

dV dV 

- — dipx = -- — dip 2 (83) 

Oipi Oip2 

on this path, one is led to the result 



1 d 4 V 

6 dff 



1 

2A2" 



d 3 V 
dpi 



(84) 



J cl 



which is compatible with, and generalizes, (pl|). 

Let a2A and h\ be the values 02 and h take as a func- 
tion of g and the other variables on the A-line. These values 
are fixed by the first one of the conditions ([Tq), Ai = 0, 
and (|8(]). To determine the coexistence surfaces p and er, 
one must again exploit the equality of the grand potentials 
of the corresponding coexisting phases [cf. (pij)]. 

Since §\ now does no longer vanish, these conditions 
are not easy to handle analytically. In general, one must 
recourse to numerical methods. Fig. ^ shows an example 
of a phase diagram obtained in this fashion for the choice 
of interaction constants given in the caption. As expected, 



a- 2 



3.5 - 



2.5 



1.5 




Fig. 5. Phase diagram of the nonsymmetric model in the Lan- 
dau approximation. The results, obtained by numerical (and 
analytical, see text) solutions of the equations of state, corre- 
spond to the following choice of parameter values: i>2 = —6, 
04 = 64 = 1/6, d,2i = —1/2, du — 0.2. Only the neighborhood 
of the critical end point is shown. 



the coexistence surface p, and hence the critical end point, 
are displaced from the h = plane. Further, p gets curved, 
and the reflection symmetry of the phase diagram with 
respect to p we had in the symmetric case is lost. 

These features can also be verified by analytic calcula- 
tions. To this end, we will restrict ourselves to the vicinity 
of the critical end point and treat the symmetry-breaking 
term tx du in a perturbative fashion. In order to determine 
the A-line, we take 62 < and g > g e , with g e = g e (dn), 
the (as yet unknown) value of g at the critical end point. 
We write 

M9) = ^x ) (9)dn + 0(d 2 n ) , (85) 



Mg) = Mg) + Ax(g) + o(d 3 n ) , (86) 



and 



a 2X (g) = \d 21 \^+(g) + 0(d 2 11 ) 



(87) 



where ip+ (g) was defined in (J4J) . The values of these quan- 
tities for du = are known from our analysis of the 
symmetric critical end point in the previous subsection. 
That the term linear in du of ^>\{g) vanishes follows from 
the fact that 4>\ and <f>\ are coupled via terms of order 
dn4>\ ~ 4>\ ~ Upon substituting these expansions 
into the condition ( [sol) and the expression (^6|) for tan d\ , 

we can solve for <j)^{g), obtaining 



Wis) 



21 



2a 4 [3 64^+(ff) 



4i 
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_ -|<W+'(g) 

2a^+{g)] ' 

where [tp + (g)] and ip+'(g) are given by © and ©, 
respectively. 

The above expansions can now be inserted into the 
equation of state (pj) to determine h\ . This gives 



Utilizing the result 
coefficient 



h x (g)^d 111 p+(g) + 0(dl 1 ) . (89) 
~"|) for ipx> one gets the second-order 



^\g) = -d\ 1 ^ + '{g)4> x 1) {g) 



1 - 



121 



(.9) 



(90) 



In order to be able to determine g e from the coexis- 
tence condition 



V[<i>\, ipy,a 2 \(g e ), hx(g e ),g e ] 
= V[4> a ,ip a ;a 2 \(g c ) 1 h\(g e ),g e 



(91) 



we must also know ip a and <p a to the appropriate order 
in G?n. The corresponding series expansions are analogous 
to the those of tpx and 4>\ given in ( pf^ ) and (|85|). The 
zeroth-order term of i/> Q is "0-(<?) = — ip+(— 5) [cf- (|38|)]. 
Hence we have 

[/» - dn l/vJ A = d n VP+{g) + ^+{-g)\ + 0(d 3 n ) (92) 
and 

[a 2 + d 21 ^ a ] x = |dai| ty+{g) + + 0(d?i) ■ (93) 

To obtain (/>„on the A-line, we insert these results into the 
expression (|66[ ) for 770, which we ex pan d in powers of its 
last argument, [h — dn ip a ]\> using (p6|). This yields 



which leads to 

$<A X = -ip+(-g) - 



\d2i\ 

d? n 
2\d 21 



0(d 3 n ) 



(94) 



i, + '(-g)+0(d 3 n ) (95) 



upon insertion into the expression ([37]) for r\- . 

Using the above results, one can determine g Q (dn) to 
second order in d\\ in a straightforward manner from the 
condition (|9l|). The result is 



9e(dn) = 



d 2 n 
2\d 21 \ 



0(d 3 n ) 



(96) 



As a check, let us set the interaction constants to the 
values b 2 = —6, a 4 = 64 = 1/6, d 2 % = —1/2, dn = 0.2 
utilized in our numerical analysis (see Fig. | ) . Our pertur- 
bative analytical results (^6|), (89), and (|7) then predict 
the critical end point to be located at 



g c ~ 0.04 , 
h e = h\(g c ) ~ 1.2 , 
a 2c = a 2 \(g c ) ~ 3.0 , 



(97) 



and to have the slopes 

a 2 \(g e ) 



d n 1>+{0*) 0.017 , 

|d 21 |V+'G?e) ^0.04. 



(98) 



These numbers are in excellent agreement with the nu- 
merical results shown in Fig. ^. 

Next, let us calculate the equilibrium densities <j>pf 
and of the ordered phases (3 and 7 in a perturba- 
tive manner. To this end, we consider a path in the space 
{{a 2 , h,g)} with g = const that is asymptotically parallel 
to the coexistence surface p and intersects the A-line at a 
point Q = (a 2 \(g), h\(g), g) close or equal to the critical 
end point. Taking Q as expansion point in (|77|), we express 
the deviations S(f> = tj> — <p\ and 5ip = ip — from the 
associated classical solution {4> c h^c\) = (0a, "0a) in terms 
of the eigendensities ip% and ip 2 of Q. 

The shifts Sa 2 = a 2 — a 2 \(g) and Sh = h—h\{g) induce 
the following changes of V(0, V0 : Nonvanishing terms 

^a»(* A )rf = ^[l + 0(<&)]<P?, (99) 



- Sa 2 ip\ip 2 sin2??A ~ dn Sa 2 (p\(p 2 



(100) 



f 2 [{<px Sa 2 - 6h)sm&x] ~ du <f 2 (101) 

are generated; the coefficient \ 2 (4>\, ip\) of the \ ip 2 2 term 
is changed by an amount 5a 2 0(d\ 1 ). The term linear in 
ipi is negligible sufficiently close to Q since we required 
the path to be asymptotically parallel to the coexistence 
surface p. 

Minimizing with respect to tp± and tp 2 yields the equi- 
librium values 



±4 



\Sa 2 \ 

2 



[l + o(da 2 )} , 



ip 2 fi = Sa 2 0(d 11 ) , 



from which we get 



±1 



' 1^2 | 
64A 



[1 + o{5a 2 )\ 



(102) 
(103) 

(104) 



and 



tpp - i)\ fa ip+'(g) 



1 , 1 \Sa 2 \ 



«4A 



\Sa 2 1 

«4A 



(105) 



Since <5a2 ~ (— t), the result displays the expected |t|^ 
singularity of (0), with the mean-field exponent /3 = 1/2. 

In order to determine the wings of a-(3 and a-7 co- 
existence of the first-order surface er, we would have to 



6 Note that Sh ~ dn. Its zeroth-order term in dn vanishes 
because the coexistence surface is given by h — in this (sym- 
metric) case. 
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substitute the classical equilibrium solutions <p p and ip p , 
with p = a, (3, 7, into V(<p,tp) to obtain the correspond- 
ing grand potentials A p [a 2l h, g) and then solve the coex- 
istence conditions 



A a {a 2 ,h,g) = Ai(a 2 ,h,g) , l = 0,j, 



(106) 



e.g., for g as a function of a 2 and h. On the triple line 
conditions (106) hold for both I = a and I = (3. 

Let us consider small deviations from the critical 
end point that are directed along the triple line. We 
parametrize these through Sa 2 = a 2 — a 2c and denote the 
value of the effective <p>\ interaction constant (Q) at the 
critical end point as hi c . Using the above results and ne- 
glecting terms of order d\i, it is not difficult to show that 
the triple-line value g a behaves as 



g a « - [1 + 0{d\ l )} 



(5a 2 y 



<a 4c y/\b 2 \/h 



(107) 



as da 2 ~ t — > —0. The second derivative of this expression 
with respect to t, taken at t = +0, gives us the expected 
jump singularity. (The corresponding value of g"(t) at t = 
+0 is zero.) 



Here ip^' are the mean values of ip in the pure phases 
selected by the choices g = ±0, respectively, while Sip 
denotes the fluctuat ing q uantity inside the brackets on the 
right-hand side of ( |l0§| ). Hence the correlation function 
becomes 

(Sip(x) Sip(x')) = w 2 S(x - x') (109) 

for \x — x'\ >£</,. 

If we let B — > in Tiity], so that ^ — > 0, and ignore 
the cubic and quartic terms in (^|), we obtain w 2 = 2\b 2 \ 
and — ±y / \b 2 \/b4. Inclusion of these terms (e.g., by 
means of perturbation theory) does not modify this large- 
scale form of the theory's correlation functions, but pro- 
duces different values of w 2 and ■ In other words, the 
parameters of this Gaussian (fixed-point) Hamiltonian are 
changed, albeit its form remains the same. In RG theory 
such a change of coordinates of a fixed-point Hamiltonian 
in parameter space is known to be caused by redundant 
operators pp . 

For the sake of clarity, let us briefly recall the concept 
of redundant operators. Consider a field theory with a 
Hamiltonian Ti.[(p], where tp could be a single-component 
field, such as <p or or a multi-component field (</? a ), 
such as ((f>, ip). Suppose H*[<p] is a fixed point under RG 
transformations. We consider operators of the form 



4 Beyond Landau theory 

In the previous section we saw that Landau theory yields 
for both the symmetric and nonsymmctric versions of our 
model phase diagrams with the correct features. We also 
verified that the results of this approximation are consis- 
tent with the predicted singularities (|l]) and (||). We now 
wish to extend the analysis beyond Landau theory. 

It will be helpful to make a number of remarks before 
we embark on details. Our strategy will be to relate the 
critical singularities on both the critical line and at the 
critical end point to those of the standard </> 4 theory. In 
other words, we will show that they can be described by 
a Hamiltonian of the form (||) of Hi [<f>]. To understand 
how this works, it is useful to consider first the theories 
described by the Hamiltonians T-ii[(j>] and TC 2 [ip] in the 
absence of any coupling between <fi and ip, and then discuss 
what happens when Ti.i 2 [cj),ip] is turned on. 



0= d d xO(x) , 0(x) = 0[(p(x),V(p(x)] , (110) 
Jn 

i.e., 0{x) is a local operator depending on tp and its deriva- 
tives, but not explicitly on x (to ensure translational in- 
variance.) Such an operator is called redundant if a local 
functional T(x, [<p]) = T[ip(x),W(p(x)} exists such that it 
can be written as 



(111) 



where C5tra[^] is the generator of transformations of TL 
induced by the change of variable 



<p(x) = <ff(x) + r(x,[<ff]). 



That is, 



(StraPI «M = l d<lx 
JO 



r(x 



5H[ip] sr(x) 



5<p(x) Sip{x) 



(112) 



(113) 



4.1 The decoupled noncritical theory and redundant 
operators 

Suppose that B > 0, b 4 > 0, g = ±0, and b 2 < b 2c , where 
b 2c is the critical value below which the symmetry ip — > — ip 
is spontaneously broken. (At the level of Landau theory, 
b 2c = 0, of course.) Then 7i 2 [ip] describes a massive field 
theory whose statistical properties on length scales large 
compared to the corresponding correlation length £^ may 
be characterized by a Gaussian probability distribution 
with Hamiltonian 



where the second term in brackets results from the Jaco- 



— H G [5iP\ = — 



tp — ip 



(±) 



(108) 



bian of the transformation (112). 

Now, con sider the noncritical Gaussian Hamiltonian 
( |l3"| ) or ( |10§| ). This is a fixed point under RG transfor- 
mations. Wegner (see Sec. III.G.2 of Re f. Jt3|) has shown 
that any operator O of the form (lit) — i.e., the inte- 
gral of any local operator — at this Gaussian fixed point 
can be expressed as a sum of a redundant operator and a 
constant: 

o = o red + c \n\ , (ii4) 

where the constant Co may be zero for some operators 
( as it must for any operator that is odd in ip). The result 
( 114 ) can be proven in a constructive manner li^ . Wc 
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content ourselves to showing how (114) reads for the cubic, 
quartic, and (V-i/>) 2 operators appearing in (|l2|); one has 



^ 3 = -0tr a [V 2 +2 5(O)]W G [^] ) (115) 
a z 

! V 4 = ^tra[^ 3 +3V5(0)]H G M-3|^|[5(0)] 2 , (116) 



in 
and 



(V^) 2 = I- ©tra[-V 2 V] WoM + 1^1 / q 2 • (H7) 
S2 Z Jq 



In a theory regularized by a cutoff A the momentum 
integration J q = J d d q/ (2n) d is restricted by |qr| < A; 

thus (5(0) = J 1 is not infinite but equal to KdA d , with 

K d = 2 1 - d n- d / 2 /r{d/2). 

It is u seful to recall the physical significance of the re- 
sult (114): The massive field theory described by 7i2[i>] 
does not exhibit critical singularities; the associated bulk 
free energy density is an analytic function of the thermo- 
dynamic fields. Thus no scaling operators other than the 
trivial one \D\ appear; the associated scaling field is the 
previously mentioned special one with RG eigenexponent 
d. 



4.2 The decoupled critical theory 

We next turn to the field theory described by the Hamil- 
tonian TLi[<p\. For given positive values of A and 014, this 
has a critical point at (aa, h) — (020 0), where a 2c = 
in the Landau approximation. The asymptotic behavior 
at this critical point can b e analyzed by means of well- 
known RG methods 



|,^fj. In order to establish our 

notation, it is necessary that we briefly recapitulate some 
of its ingredients. 

Let us introduce the corresponding bulk free energy 
(per volume and k^T) 



f 1 . h (a 2 ,a 4 ,h) = - lim (r^ln (v^e~ H ^ 



(118) 



as well as the generating functional Z\ [J, I] and Gi [J, I] 
of correlation functions and cumulants, respectively, via 

Z 1 [J,I] = e g ^ = /e««+K^ 2 )\ , (119) 

\ / 'Hi 

where the subscript Hi serves to indicate that the average 



,7il ~ /Xtye-«i[ 



(120) 



refers to the decoupled system with Hamiltonian 
(rather than to the full one with Hamiltonian TL[4>, ?/>]), 
and 

(I,cj) 2 ) = [ d d xl(x)cj) 2 (x) (121) 
Jo 

is a convenient shorthand. 



To absorb the ultraviolet (uv) divergences of the theory 
ind = 4 — e<4 dimensions, we use reparametrizations of 
the form 



<P = [z<f,(u)] 1/2 4> mii , 

a 2 = ft 2 Z T (u) t + a 2jC , 
d4 = K e Z u (u) u , 

h = K ( d+ 2)/2 [Z<ji{u)] -l/2 hie 



(122) 

(123) 
(124) 

(125) 



Here k is an arbitrary momentum scale. The renormal- 
ization factors are understood to be fixed my means of a 
r ('mass') independent renormalization scheme; for con- 
creteness, we shall assume that dimensional regularization 
is employed and that they are determined by minimal sub- 
traction of poles in e. In a cutoff-regularized theory, the 
shift a2, c of the critical point from its zero-loop value 
diverges as A 2 as A — > 00. 



In addition to the reparametrizations (|122|)-(125), ad- 



ditive counterterms are required. The bulk free energy /i^ 
and its first and second derivatives with respect to a 2 have 
primitive divergences ~ A 4 , ~ A 2 and <~ ln^l at d = 4, 
respectively. To cancel these, we make subtractions at the 



normalization point (NP) a 2 



,NP 



h = 0, with 



< P =a 2 , c + K 2 Z T r NP 



r NP = l 



(126) 



Accordingly we introduce the renormalized bulk free en- 
ergy 



A re b ( r ' u > = A,b( a 2, a 4 , h) - / ljb (a 4 ) 
-(«2-«f)e NP W 

+^(« 2 -«n 2 ^ NP M, (127) 

where the renormalization functions e NP and C NP are 
fixed through the normalization conditions 



/r b n (i, u ,o) = o, 



3/r b n 



and 



dr 

-^-(l,u,0) = 



(1,«,0) = 



(128) 



(129) 



(130) 



For the generating functional of renormalized cumulants 
we have (cf. Sec. 12.1.1 of Ref. @) 



GT n [J,I]=Gi 



^6 ^ J, Z T I 



e NP Z T I 



\ C ^[Z T I-(a 2 -ar)} 2 ). (131) 



Thus the renormalized cumulants 



r (N,M) 
l.rcn 



5 N+M GT n [J,i] 



5J{ Xl ) . . . 5J{x N ) SI{Xi) . . . SI(X M ) 



.1 = 
1 = 



(132) 
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are related to their analogously defined bare counterparts 



n (N,M) 
"-^l.rcn 



ZI N ' 2 Z. 



■M ^q(N,M) 

5%l [e NP -C NP K 2 Z T (r-l) 



(133) 



where Xi 2 = X\ — X 2 . 

Upon varying k, we can derive RG equations. Let us 
introduce the beta function 



P u (u) = nd K \ Q u , 

the exponent functions 

r)<j,(u) = Kd K \ n InZ^ , 

■q T {u) = k8 k \ q \nZ T , 
and the differential operator 

V K = K d K +/3 a d u -(2+r] r )Td T - d+2 ~^ 



(134) 

(135) 
(136) 

hrendh r<!a , (137) 



where <9 K |o means a derivative at fixed values of the bare 
parameters a 2 , 04, and h. Then the RG equations can be 
written as 



NP 



(138) 



and 



N 

V K + —T ltj) -Mr lT 



r (N,M) 



- \PtJX ,b) Vl ( ) 



,.2M 



M.2 



(139) 



The RG equations ( 



( 13q ) and (139) can be exploited in 
a standard fashion to derive the familiar scaling forms of 
the (singular part of the) free energy and of the cumulants 

G\ N '*! } - Details can be found, for instance, in Refs. |48|, 
|49| , l50t or Q and will not be repeated here. In the 



cases 

of the free energy f[ e £, the energy density e rcn = G^^, 

and (specific heat) the RG equations are inhomo- 

geneous. Explicit solutions to such equations in terms of 
RG trajectory integrals are given in Ref. (52) (and else- 
where). Specifically for the renormalized free energy /[™, 
the solution reads 

fl^(r, h len , U, k) = (n£) d fi$[f(£), h(£), u(i), 1] 

+ /^^)- ( 14 °) 



with 

m 



[E T (u(l),u)Y 



(141) 



where u(£), t(£), and h{£) are standard running variables 
[cf. Eqs. (3.79a), (3.79.b), and (3.79d) of Ref. [|l|]. Utiliz- 
ing the notational conventions of this reference, we have 



'-(£) = e- 1 ^ E t [u{£),u]t r> e- x ^E*( u ) 



and 



h(£) = f A l v E h [u(£), u] h ~/~ A,v K( u ) h ™ 



Here 



and 



E T [u, u] = exp ■ 



Eh[u,u] = exp ■ 



r u{i) 



dx 



Vr - Vr(x) 
Pu(x) 



r u{l) 



dx 



(142) 



(143) 



(144) 



(145) 



An asterisk is used to mark values at the infrared-stable 
fixed point u* of the beta function (3 U . Further, E* h (u) 
means E T .h(u* , u), and the exponents are given by v — 
(2+r]*)- 1 and Ajv = (d+2-r$)/2, respectively. 

The consequences of these equations for the free energy 
/[™ can be cast in the usual form: It is a sum of a regular 
and a singular part, 



rroii _ rreg . /-sing 
./l,b — ./lb Jl,b 



(146) 



The latter, a solution of the homogeneous RG equation of 
/{ ", near criticality takes the scaling form 



f^(r,h ien ,u)^\fr a Y ± (h\f\- A ). 



(147) 



Here f ss E*(u)t and h sa E^(u) h lcn are the two rele- 
vant scaling fields. The scaling function Y±(y) is universal 
(and hence independent of u); its two branches ± satisfy 
matching conditions as y — ► ±00 (see, e.g., Ref. Q). 



4.3 The coupled theory 

We are now ready to discuss the coupled theory, with 
7i\ 2 ^ 0. We begin with the symmetric case, setting 
h = d±i = and A = 1. To simplify the subsequent 
analysis, we also set the interactions constants e2i and 
/21 to zero, but will comme nt on their effects later [see 
the paragraph following Eq. (163)]. 



4.3.1 Symmetric case 

As reference densities about which we expand we take the 
classical solutions <j) c \ — and if) c \ at a point (0,2, h,g) = 
(ajvv (<?ci)j 0,<7 re f) with g rc { > g e on the Landau-theory A- 
line. Writing 



tp = tpcl + tp , 



(148) 
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we have 



with 



H[<p,il>] =H 2 bPci]+H[cf>,i>]+H'l({>,i>} 



(149) 



with 



H[<t>J] 



2 V -r 2 v V 



(150) 



and 
where 



a 2 = a 2 

6 2 = b 2 

63 = 36 4 V'cl 



C?21 V'cl 

I 

'l ' 



36 4 v> 2 



(151) 

(152) 
(153) 
(154) 



Suppose first that 7i' is neglected. Since Ti is quadratic 
in "0, having the same form as the Hamiltonian of the 
stochastic dynamic models C, D, and E ]44}| , the field ^ 
can then be integrated out exactly to obtain an effective 
(j) A model. Let us define quite generally both for Ti' = 
as well as for Ti' ^ 0, an effective Hamiltonian through 



-H. 



off LP] = 



Vipe 



-H 2 Wci]-ii[(l>,#]-H'[4>,i>] 



For TV — we have 



Heft [ 



1 . &2 
= X m — 

2 2tt 



7Yi[^;a 2 , 



(155) 



(156) 



r 3 = ^- 



2 



2b\ 



and 



2b\ 



d 2 i ■ 



+ 8 6 2 <5(0) 



(160) 



4k 



12b 2 S(0) 7 



Ab\ 



d 2 



8 6| 



[20 6 2 c5(O) +dj 



(161) 



The meaning of these results is the following. To first 
order in the coupling constants 63 and 6 4 , the terms 
oc J Q ip 3 and oc J n tp 4 can be transformed away at the ex- 
pense of (i) additional irrelevant interactions (oc cf> n , with 
n = 6,8), (ii) contributions oc \Q\, and (iii) a change of 
the coefficients of the 4> 2 and <fi 4 terms. 

The irrelevant interactions may be dropped, (iii) can 
be absorbed through a redefinition of the temperature- like 
scaling field and the irrelevant scaling field ~ u—u*, which 
hence become dependent on 63 and 64 (as well as on ei 2 i 
and 6 2 ). (ii) means contributions oc 63 and oc 64 to the con- 
stant part of the Hamiltonian, which we write as /iq \ {2\, 
where /iq is Wegner's special scaling field (cf. Sec. III.G.2 
of Ref. |Q). Since 6 2 depends via i/> c i on the nonordering 
field g, these contributions to /zrj have g- dependent parts. 

We are now ready to understand the origin and nature 
of the discontinuity eigenperturbation indicated in Fig. 0. 
Let fj,Q (g) = [aaA (<?) 1 <?] be the value of (mq at a point 
(a 2 x(g),h=0,g) on the A-lin^, and /i{, o) = ^ (sO, the 
corresponding critical-end-point value. Consider a varia- 
tion g c — > g c + 5g, a 2o — ► a 2c + Sa 2 along the A-line. This 
induces a change 



where Ti 1 [</>; a 2 , a% s ] is the h = variant of the Hamilto- 
nian (|^), with a 2 and 04 replaced by expression (152) for 
a 2 and 



04 



4i_ 

2 6, 



(157) 



respectively. Thus, if TV can be set to zero, the critical 
behavior of the Hamiltonian (14E) reduces indeed to that 
of a 4> A model. This result is, of course, not new: it has 
been used in the construction and analyses of the dynamic 
models C, D, and E for a long time (Cf. also Ref. 

Next, consider what happens when TV is included. A 
straightforward calculation shows that the cubic and quar- 
tic terms of TV can be rewritten as 



1 a 
and 



4i_ 

8 61 



3d 2 iJ(0) , 2 
-<P 

(158) 



2 61 



^ = e5 tra [r 4 ]H[0,v]- 

n Jn 



4i 



16 6| 



+ 



3^ 5(0) 



2 61 



(159) 



6^=^\g)-^ 



of the special scaling field //o, with 5/i^ ~ T c (g) 

for small Sg. The shift of the Hamiltonian Sfi^ \£2\ is an 
eigenperturbation with eigenexponent y — d. Obviously 
this is the one we were looking for and whose eigendirec- 
tion is shown in Fig. |[ 

We still have to consider the (V-0) 2 term of TV. Power 
counting tells us that it is irrelevant. One also verifies that 
it can be rewritten as 



(162) 



(VVO 2 = <S tra [-b 2 - 1 Ai>] fi[<]>, $ 



d 2 \ 



(163) 



The last term in the second line contributes again to the 
constant part /zo |!?| of the Hamiltonian. The first one, 
proportional to J_ <fi 2 Atp, is irrelevant according to power 

7 We suppress the dependence of fio and a,2\ on the other 
variables, 62, 64, etc. Instead of the variables ai, . . .one could, 



of course, also use 02, 
of the former. 



Yet we prefer to express /io in terms 
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counting, w hich conforms with our expectation that the 
interaction (163) is irrelevant (apart from contributions to 
the constant part of the Hamiltonian) and can be dropped. 
Finally, power counting also indicates that the terms in- 
volving e2i and /21 are irrelevant, so that they may be 
omitted as well. 

Having iden tified the sought discontinuity eigenper- 
turbation (162) and knowing that it does not contribute 
to the critical singularities on the A-line or at the end 
point, we can now focus directly on these, droppin g th e 
ip 3 and tp 4 inter action s as well as the (V^) 2 term (163). 
From Eqs. ( |155| ) and ( |l56| ) we know already that the re- 
sulting truncated Hamiltonian Ti. reduces to an effective 
<p Hamiltonian upon integrating out ip. Moreover, the 
Hamiltonian Ti and its renormalization are well known 
from studies of the stochastic dynamic models C, D, and 
E (see, e.g., Refs. Q, ||, or Sec. 35.4 of @). Adding 
a source term —JqIiP to Ti, one can compute the func- 
tional integral §Vip exp{— Ti [(f), ip] + J n Iip} to relate the 
cumulants (Iljli 4>( x j) YlfLi ^(Xi))^" 1 pertaining to Ti 
to cumulants of a standard <p 4 theory described by the 



Hamiltonian (156). One has 



^(X))^ = -^-(p 2 (X)) Hc!{ (164) 
zo 2 



and 



b 2 4 b% 

(165) 

All other correlation functions of the fields ip and <p are 
identical to their analogs of the effective (p 4 theory one 
obtains through the replacement 



2b 2 



(X) 



(166) 



Thus the critical behavior of the above cumulants 
(IljLi 4>( x j) Ili^i ^{Xi))^" 1 may be inferred directly 
from the known one of the corresponding cumulants of the 
cp 4 theory. Alternatively, one could derive RG equations for 
the former, using the reparametrizations Z u (u) u, 



tp = Z l J 2 \u) 4> mn , ip = b. 
n 2 Z t (u)t + a c 2 s c , a 2 



I-1/2 7 i/2 

J 1> 



(u,j)ip len , a' 



k 2 Z^{u,~t) Z T (u) f + a 2jC , and 
62 k"/ 2 z]/ 2 (u,j) Z T (u)j. Here Z u , Z^, and Z T 



are t he renormalization factors introduced in Eqs. (122)— 
( p5| ), while Z$ is of the form Z$(u, 7) = 1/[1 - 7 2 f[u)], 
where f(u) is known to be given by t he L aurent part of 
the additive counterterm oc I 2 in Eq. ( |l3l[ ) ||||5§]. 

From Eq. (131) we see that the leading singularity of 
(tp) (and hence of (ip)) on a path with t ~ T—T c (g) agrees 
with that of the energy density (<p 2 ). Since on the liquid 
branch of the liquid-gas coexistence boundary the scaling 
field t varies asymptotically linearly in t ~ T—T e as the 
critical end point is approached, the leading singularity 
of the density g tot — (ip) is indeed of the form V± |£| 1_Q , 
where is the usual universal specific- heat ampli- 

tude ratio. 



4.3.2 Nonsymmetric case 



Considering points (a 2 ,h,g) on the fry (disordered) side 
of the coexistence surface p, we expand about the clas- 
sical solutions cpp and ipp~- Since we ultimately wish to 
approach the A-line or critical end point, we may restrict 
ourselves to points in their immediate vicinity. Introduc- 
ing the shifted densities ip and ip via 



p> = p> c \ + p> , tp = ip c i + ip, (167) 

we express the Hamiltonian defined through Eqs. (^) and 
(H) - ® m terms of these. The resulting quadratic ('Gaus- 
sian') part can be written as 



Vi?-r7 cl A 



ip 



with 



*7ci 



and v£ } 



A + /21 ip c \ en + e 2 i 4> c \ 
en + e 2 i <p c \ B 

;( 2 



(168) 



(169) 



ipc\), where V is the matrix intro- 



V (2)( 

duced in Eq. 

The matrix ?7 c j must be positive definite. [Otherwise, 
the classical homogeneous states (4> c h V'ci) are not guaran- 
teed to be stable.] It plays the role of a metric. Hence we 
can diagonalize Tic by a similarity transformation 



iP 



= U 



(170) 



satisfying the diagonalization condition ( [74] ) together with 
the orthonormality relation 



U l Vcl U = l 



(171) 



to obtain 

1 



n, ; - - jj^^i) 2 + + (V^ 2 ) 2 + A 2 f 2 2 ] . (172) 



By analogy with Eqs. ( |149| )-( |151[ ), we split the resulting 
total Hamiltonian (||) as 

H[<l>, ip] = H[4> ch ip d ] + ft[<pi,<pa] + H'[ip u <p 2 ] (173) 
with 



Ai 



A 2 2 

Y <P2 



and 

H'[<pi,<p2] = 



2 
V2.1 



fi <P2 



(174) 



~fv>2 
J=3 J 



+ Ul,2 ¥2 + v 3,l <Pl f2 + V 1>3 if 1 lfi 2 
'±j> 2 Al/j 



V2.2 <p\ H>\ 



(175) 
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We have dropped a contribution U3.0 <^ 3 to the inte- 
grand of 7i'; its coupling constant 1)3^ vanishes at the 
Landau-theory location of the A-line by condition (|79|), 
and we know already that it can be transformed away by 
means of a shift ip% — > ipi +^1; i.e., aside from a change 
of the interaction constants of the Hamiltonian Ti. and its 
constant part, it corresponds to a sum of redundant and 
irrelevant operators. 

The Hamiltonian H has the familiar model-C form 



( |l50| ). The contributions in the f irst line of Eq. (175) are 
the analogs o f the Hamiltonian ( |l5l| ); from our consider- 
ations in Sec. 4.3.1 we know their effect: Apart from con- 
tributing to the special scaling field (and shifting the A-line 
and critical end point), they involve redundant and irrele- 
vant operators, and hence may be dropped when analyzing 
the asymptotic criti cal b ehavior. The same is true of the 
remaining terms of (151). The easiest way to arrive at this 
conclusion is via power counting, utilizing the canonical 
scale dependences ipi ~ k 1 " 6 / 2 and (f2 ~ n 2 ~ e / 2 . More- 
over, the contributions involving the interaction constants 
ui,2 , ■ • ■ ,V2,2 could be analyzed along the same lines as 
the Tp 3 and ip nonlincarities in the previous subsection by 
choosing appropriate functionals T[y>2] to eliminate them 
via a change of variable ip% — * ^2 + ^l^]- The contribu- 
tions proportional to e2i and /21 correspond to a variety 
of derivative terms ipfAipi, <p 2 A</?2, ■ ■ ■ , (Vy^) 2 ^, all of 
which are irrelevant according to power counting. 

Upon ignoring H' altogether, we get back to the 
Ha milton ian H[(pi, ^2], which is equivalent to the one in 
Eq. (150). The implications for the asymptotic critical be- 
havior that occurs when the A-line or the critical end point 
are approached is obvious: ipi, the order parameter, and 
the secondary density tp2 behave asymptotically as their 
respective analogs <f> and ip m the symmetric case. 

The density ip has components along both ipi and if2- 
Consequently its average gtot = (ip) has singularities of 
the form U% \tf (with U\ = 0) and V± in the limit 

t = (T — T c )/T c -> ±0. This holds, in particular, for ap- 
proaches along the coexistence boundary g a {t). Thus the 
leading singularity of the liquid density g t ot is ~ |i| ^ 
the critical end point is approached along the triple line, 
and ~ lil 1 " if it is approached along the fij-a section 
of the liquid-gas coexistence line. These findings are in 
conformity with Eq. (0). 



4.3.3 The singularity of the first-order line g a {t) 

A crucial assumption on which the phenomenological 
derivation || of the |i| 2 ~ Q singularity of the first-order 
line g a (t) is based is the asymptotic behavior of the grand 
potentials A p of the phases p = f3, 7, and fij at the criti- 
cal end point. This is hypothesized to be of the same form 
as at a conventional critic al po int; i.e., it can be decom- 
posed by analogy with Eq. (146) into a regular background 



contribution and a nonanalytic part, A s J, ng , which has the 



asymptotic scaling form ( |147[ ). 

Keeping in mind our identification (162) of the dis- 
continuity eigenperturbation, one realizes that our above 



field-theoretic considerations are in full accordance with 
these assumptions, corroborating them. Utilizing them 
one can proceed |39f | just as in the derivation [g| based 
on the phenomenological scaling theory to determine the 
coexistence singularity of the first-order boundary g a {f) 
from the coexistence conditions (p6|). Since f ~ t for small 
deviations from the critical end point directed along <?o-(i), 
this line inherits the |i| 2 ~ a singularity of the grand poten- 
tials Afjj and Ap y. Thus Eq. (|l|) follows. 



5 Summary and conclusion 

To put things in perspective, it will be helpful to reca- 
pitulate the main steps of our analysis and summarize its 
principal results. 

(i) In order to study the critical behavior at symmetric 
and nonsymmetric critical end points as well as the asso- 
ciated coexistence singularities, phenomenological contin- 
uum models with a Hamiltonian depending on two fluc- 
tuating densities <p and ip were introduced. These models 
may be viewed as continuum variants of the BEG model 
[p3| . They generalize similar ones previously considered in 
the literature (cf. Ref. fl3fl), to which they reduce when 
some interaction constants vanish. We justified them via 
phenomenological arguments, but were also able to de- 
rive them directly from the BEG model, utilizing a gener- 
alized Kac-Hubbard-Stratonovich transformation to map 
the latter exactly on a lattice field theory and a subsequent 
continuum approximation. 

(ii) We then employed Landau theory to analyze the 
models for the cases in which the Hamiltonian obeys or 
breaks the symmetry <p — > — <p. Provided the interactions 
constants were chosen in the appropriate ranges, phase di- 
agrams with either a symmetric or nonsymmetric critical 
end point could be obtained (cf. Figs. g-g). Furthermore, 
we corroborated that the |i| 2_Q singularity ([!]) of the first- 
order line g a (t) and the |i| 1_Q and \t\" coexistence singu- 
larities (Q) of (ip) take on the expected Landau-theory 
forms, namely, discontinuities in the second derivative of 
g a (t) and in the first derivative of (-0), or a singularity 
~ \t\ x ^ 2 , respectively. 

(iii) Building on these results, we then studied the crit- 
ical behavior on the A-line and at the critical end point. 
One important observation was the following. In order to 
obtain phase diagrams with a critical end point (in Lan- 
dau theory and beyond), it was essential to include non- 
linearities such as a ip 4 term that are irrelevant according 
to power counting. The presence of this term guarantees 
that a stable a (gas) phase exists, even in the absence 
of coupling to the (primary) order-parameter density <p- 
However, once it has been verified that the phase diagram 
has the correct topology (with a critical line A, critical 
end point, and first-order boundary g a (t)), one can focus 
on the more specific problem of the asymptotic critical 
behavior that occurs when either a point on the A-line 
or else the critical end point itself is approached along a 
given thermodynamic path. Then approximations that are 
tailored for this particular purpose become viable. 
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Upon expanding about the classical equilibrium values 
of the densities <f> and "0, and expressing the Hamiltonian 
in terms of the deviations <f> and ip from these, the Hamil- 
tonian was found to involve also nonlincarities that are ir- 
relevant according to power counting, such as the ip 3 and 
?/> 4 terms in the symmetric case. We were able to show, 
within the context of perturbation theory, what their ef- 
fects are: Aside from (a) inducing a change of the usual two 
scaling fields (the temperature-like scaling field r and the 
magnetic-field like scaling field), they (b) correspond to 
irrelevant operators (giving rise to corrections to scaling) 
and (c) contribute to the constant part of the fixed-point 
Hamiltonian, i.e., to Wegner's Q special scaling field /xo- 

(iv) These findings enabled us to verify the existence of 
a discontinuity eigenexponent d at the fixed point repre- 
senting the critical end point, and to understand its origin: 
The eigenperturbation with which it is associated corre- 
sponds to the change dfio of the special scaling field no 
that occurs when one moves away from the critical end 
point such that the theory remains critical. 

(v) In conjunction with the findings mentioned at the 
end of (iii), (iv) shows that the critical behavior on the 
A- line is the same as at the critical end point, inasmuch 
as the values of the critical exponents and other universal 
quantities are concerned. 

We emphasize that our conclusions mentioned in (iv) 
and (v) are in conformity with the position-space RG re- 
sults of Ref. [§2|,|2^] and the RG picture that has emerged 
from them (cf. Fig. ||). Our results show how this RG pic- 
ture translates into field theory; they verify the existence 
of the discontinuity eigenexponent, clarify its origin, and 
reveal how its presence can be reconciled with the antic- 
ipated equality of the critical exponents associated with 
the critical line and the critical end point, respectively. 

(vi) Once we had established the results (iv) and (v) , 
we could derive the coexistence singularities (|l]) and (^J) in 
a fairly straightforward fashion. The singularities in 
Eq. (Q) are due to the coupling of the secondary density, 
ip, to the energy density. In the case of a symmetric crit- 
ical end point, this is the leading temperature singularity 
of (tp) on the liquid branch of the liquid-gas coexistence 
curve. In the case of a nonsymmetric critical end point, 
the densities <fi and ip mix. Thus the order parameter be- 
comes a linear combination, ipx, of both, so that if) couples 
directly to ip\. This implies that the leading coexistence 
singularity of (i/j) on the three-phase-coexistence side of 
the transition becomes ~ \t\°. The |t| 2 ~" of the first-order 
boundary go-(t) follows in a way known from its original 
phenomenological derivation |5| from the usual free-energy 
singularity of the j3j, j3, and 7 phases at the critical end 
point. 

(vii) A challenging issue raised by Fisher and Barbosa 
in their critical assessment || of the theory of critical end 
points, is the potential occurrence of essential singularities 
at first-order boundaries like a. Since our approach relied 
on perturbative RG arguments, essential singularities are 
beyond its scope. 

The following should also be clear: We cannot, of 
course, rule out that for special models the critical be- 



havior at the critical end point might differ from that on 
the critical line. However, this should not happen in the 
generic cases we were concerned with here. 

(viii) The present work, which was focused exclusively 
on bulk properties, provides a basis for investigating the 
problem of critical adsorption at a noncritical a-/?7 inter- 
face. This issue will be taken up in a separate publication. 

(ix) We close with a comment on possible implications 
of the recent work by Fisher et al |55|,[56| that appeared 
after completion of the present investigation. Upon reana- 
lyzing two-phase heat-capacity data of the one-component 
fluid propane (C3H8), these authors concluded that both 
contributions to the specific heat 

C{> ot = VT(d 2 P /dT 2 ) v - NT{d 2 fi/dT 2 ) v , (176) 

the pressure derivative d 2 p/dT 2 = p" p and the chemical 
potential derivative d 2 /j,/dT 2 = /x" p diverge as \t\~ a when 
t = (T — T c )/T c — ► —0 on the vapor-pressure curve. This 
behavior, which Yang and Yang Js?} originally suggested 
to be the most likely one for real gases, differs from that 
of simple lattice gas models, for which /i vp is found to be 
analytic through T c so that /z" p (T c — ) must remain finite. 
Building on Rchr and Mermin's earlier work ^8|, Fisher 
and Orkoulas |^| concluded that the standard scaling the- 
ory for fluid criticality must be revised inasmuch as the 
pressure difference p — p c should, in general, also mix into 
the scaling fields, in addition to the chemical potential 
and temperature differences /i — [i c and t, respectively. As 
a remarkable consequence of this 'pressure mixing' they 
found that the arithmetic mean (p\ + p g )/2 of the gas and 
liquid densities pi jg at coexistence should have a temper- 
ature singularity ~ \t\ 2 ^ as t — ► —0. Such a singularity 
would dominate over the usual energy-density singular- 
ity ~ |i 1 1 ~ a that is usually given in textbooks |I| as the 
leading one of this quantity. 

If this pressure mixing must indeed be taken into ac- 
count in the analysis of critical behavior of simple fluids, 
one anticipates it to play a role also in the case of criti- 
cal behavior at nonsymmetric critical end points. Clearly, 
the leading singularity ~ \t\P of the coexistence density 
(||) must prevail, but one expects pressure mixing to pro- 
duce a subleading one of the form |i| 2/3 , which would be 
stronger than the subleading singularity ~ |i| included 
in Eq. §. 
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for taking the time to carefully read the manuscript. This work 
has been supported through the Deutsche Forschungsgemein- 
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A Mapping of the BEG model onto a field 
theory 

We start from the Hamiltonian ( O ) of the BEG model. 
To set up the Hubbard-Stratonovich transformation, we 
introduce the row vectors 



cr =(..., Si Sj, ...) , (177) 
F =(..., Hi = H , . . . , Dj = D , . . .) , (178) 



the transposed column vector <r T , and the matrix 

' J + Mq 1 L 
L K + Mn 1 



M = ( ; T - . r , I • (179) 

Here J = (Jy), if = (ify), and i = (Ly) are the ma- 
trices of interaction constants, i.e., Jy takes the value J 
or vanishes, depending on whether i and j are nearest 
neighbors or not, and likewise for L and K . The parame- 
ter Mo has been chosen such that the matrix M is positive 
definite (cf. [Q and ]59|), so that the inverse of 7W exists. 

With these definitions the Hamiltonian ( jl4| ) can we 
written as 

Wbeg = -^o--(M-M )-<7 T -F-<7 T . (180) 

Proceeding in a standard fashion by making a Gaussian 
transformation, we can perform the trace over the spin 
variables in the partition function (fq). This gives 



Zbeg = c /o 



n/ *i* 



-Witt I^.i/-] 



with 



and 



/o = -Trln[27r(M-M )] 



(181) 



(182) 



(183) 



where 



ln[l + 2c 



iPi+D-Mo 



COS 



h(& + fl-)](,184) 



while Py, Q^, 



and are the elements of the matrices 



and 



n 



[J + M ] L K + M - L (J + M ) 



1 

Mo 



L LJ + KL 

Mo + M^ 



(187) 



respectively. Thus the BEG model (Q) is exact ly e quiva- 
lent to a lattice-field theory with Hamiltonian ( ]183] ). 

Just as in the simpler Ising case [^6,59 , the bonds 
of this Hamiltonian (i.e., the off-diagonal elements —Vij, 
— Qij, and —TZij) extend beyond nearest neighbors (NN). 
From the expansions given in the second line of (|18q) - 
(187) one sees that the NN couplings are given by J '/Mq, 
K/Mq, and L/M§, respectively, up to corrections that are 
smaller by a factor Mq -1 . Likewise, next-nearest-neighbor 
couplings are down by this factor, compared to the NN 
bonds. 

In our interpretation of the BEG model as a model 
for binary fluids, interchanging A and B particles corre- 
sponds to the transformation Si — > — Si for all sites i. 
Consequently, Hbeg is AB-symmetric if H = L = 0. In 
this case, 7iift evidently is even in <fi, as it must. 

So far our reformulation of the BEG model has 
been exact. We now make a continuum approxima- 
tion. Introducing smoothly interpolating fields (<p v (x)) = 
(<j>(x), tp(x)), we write <$? = <j> v (ai), where a is the lattice 
constant, and use the expansion 



<t>i ~ <t?j = ry-V&, [a (i + j)/2} + 



(188) 



for small = a(i — j). (To simplify our analysis, we 
take the sites i to be those of a simple d dimensional cubic 
lattice with periodic boundary conditions.) 

In this manner we obtain from Tiift the Hamiltonian 



H[</>,^] = Jd d xl^<j> v (V A VV ,W + T vv ,^ 4>v 



a w 



with 



and 



(A 



-l 



(T vv ' ) 



2da d f-^ V 

Ttoj Qoj 

j^o 



Oj 



3^0 



Oj 



(189) 



(190) 



(191) 



Q = 



J + M -L(K + M ) l L 



1 

Mo 



1 

Mo 



J 

Mo 



M 2 



K + M - L ( J + M ) x £ 



if 

Mo 



if 2 



M 2 



(185) 



(186) 



Here ^ acts to the left, while V acts as usual to the 
right. Derivatives of higher than second order have been 
dropped. 

For the BEG model (|l4|) with NN interactions and 
positive values of the NN bonds J, K, and L, the matrix 
(A vv r) is positive definite provided J + K > 2L. To see 
this note that 



1 d ^j-^-j 
a d d dq 2 



(192) 
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where 



References 



M(q) = 



Iff 



is the Fourier transform of M. But 

d \ \ dM(qY 

Mfo) = -M(g) ■ -J^- 



and 



d? M ^ 



\q=0 



M(q) 



J L 
L K 



(193) 



(194) 



(195) 



The latter is positive definite becaus e of our assumption 
J + K > 2L. In conjunction with ( |l94|) it follows that 
dM(q)/dq 2 and hence (A vv i) are also positive definite. 
Thus (A vv >) provides a metric, and the quadratic form 



T 



_ d d 2 w((f>,ij) 



Otp V 0<p V ' = 1> = Oy 

can be diagonalized by a similarity transformation 

<p v (x) = U vv > (j) v '(x) , (197) 

where (U vv i) = U is orthogonal with respect to the metric 
(A vv r): 



U' 



(A vv ,)-U=l. (198) 

For vanishing L, (A vv >) is diagonal, but for L ^ it is 
not and the quadratic part of the Hamiltonian (189) must 
be diagonalized by such a linear transformation. Thus even 
at this Gaussian level of the theory, the fields <f> and i\) 
'mix', i. e., the order parameter becomes a linear combina- 
tion of 4> and tjj, a feature that is expected quite generally 
in the nonsymmetric case L ^ 0, H ^ 0. 

To compute the interaction constants A vv >, we expand 
in L, keeping only contributions up to first order in L. This 
yields 



A = 



1 y 

2da d M Q 4r< 

3^0 



[j(M + jy 



B = 



1 



2da d M 



^[KiMo + K]- 1 } 



r 2 



and 



0(L 2 ), 

(199) 
KO(L 2 ), 

(200) 



1 + 0(L 2 



en = 



3^0 



2da d M 

+ {M + K)- 1 L(M Q + J)- 1 ] 



J)- 1 L{M + Ky 



-n r 2 



(201) 



<j> v (x)(f) v >{x) (196) 8 
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Finally, we expand the function uu((j),ijj) in Eq. ( |189| 
into powers of <f> and ip, eliminate the resulting 4> 3 an> 
terms through the shifts ( |l0| ) and (^), and drop terms of 
sufficiently high order. The resulting continuum Hamilto- 3^ 
nian then takes the form specified by Eqs. (|) and (|)~(|), 
except that the coupling constants e2i and /21 vanish in 
the used approximation. However, such interaction terms 
will be generated when one coarse-grains to a larger scale 
by integrating out the corresponding short wave-length 32. 
degrees of freedom. 
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